function [ Z ] = getKernelGradient( p, q )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

global poissonRatio shearModulus;

r = p - q;

rn = norm(r);

rg = r / rn;

M = [rg(1)*rg(1), rg(1)*rg(2); rg(2)*rg(1), rg(2)*rg(2)];
I = diag(ones(2, 1));

DrM = getDrM(r, M);
DrI = getDrI(r, rn, I);

Z = -1 * ((3 - 4 * poissonRatio) * DrI - DrM) / (8 * (1 - poissonRatio) * shearModulus * pi);
% Z = DrI;

end

